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1. Introduction 

Thermal resolution test targets are used to char- 
acterize the performance of infrared imaging sys- 
tems [1]. Low-contrast test targets are used to 
measure the minimum resolvable temperature dif- 
ference [2] that can be distinguished by these sys- 
tems. 

A commonly used thermal resolution test target 
of conventional design is shown in figure 1. The 
front plate is blackened so that it is an efficient 
emitter of infrared radiation. The back plate, which 
is also blackened, is visible through the slots in the 
front plate. The temperature of the front plate is 
allowed to float with the ambient, and the tempera- 
ture difference between the back and front plates is 
sensed by a thermocouple. An electrical current 
flowing through a heater covering the rear surface 
of the back plate is controlled to maintain the tem- 
perature difference between the back and front 
plates at a user-adjustable value. 

The simple design of figure 1 works very well 
for temperature differences of the order of 10 K 



and has been used successfully for much smaller 
temperature differences. Unfortunately, an air gap 
exists between the front and back plate causing a 
lack of tight coupling between the temperature dis- 
tribution over the back plate and that of the front 
plate. As a result, it is possible to have small lateral 
temperature gradients across the back plate that are 
not correlated with similar gradients in the front 
plate. The temperature difference at any slot is the 
sum of the intentional temperature difference set by 
the user and the difference between the uncontrol- 
lable temperature gradients. This creates a problem 
at low contrasts. 

As the temperature difference between the two 
plates is reduced, these uncorrelated gradients be- 
gin to interfere with the pattern generated by the 
slots in the front plate. For very small temperature 
differences between the plates, the uncorrelated 
lateral thermal gradients can result in the slots 
appearing hotter than the background on one 
side of the target, while appearing cooler than the 
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background on the other side. Consequently, it is 
not possible to know the actual value of the ther- 
mal contrast at any slot in such a situation. 

During the last 20 years, the performance of 
thermal imagers has improved so much that it will 
soon be necessary to test systems with thermal con- 
trasts of the order of 1 mK. It is not at all clear that 
the conventional design can meet this need. The 
purpose of this paper is to describe a new design 
specifically tailored for low thermal contrast. The 
principal advantage of this new design is that the 
lateral temperature gradients in the hotter surface 
are strongly correlated with those in the cooler 
surface as a result of the way the thermal contrast 
is generated. 




Slotted 
front plate 



Temperature 
controlled 
back plate 



Figure 1. A conventional thermal resolution test target. The 
temperature of the back plate is controlled to a constant temper- 
ature difference relative to the front plate. The temperature of 
the latter is allowed to float. 

A cross section of the new design is illustrated in 
figure 2. The target consists of an insulating mate- 
rial whose bottom surface is in good thermal con- 
tact with a heat sink, and whose other surfaces are 
well insulated. Power is applied with a uniform 



density over the top surface of the material. The 
heat is confined to flow in the vertical direction by 
the boundary conditions and the geometry of the 
device. In this case, the temperature profile over 
the top surface will follow its contour as illustrated 
in figure 2 and discussed in more detail in the next 
section of this paper. 
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Figure 2. Cross section of a target of the new design illustrating 
the principle for achieving low thermal contrast with minimum 
influence from parasitic temperature gradients. Power is dissi- 
pated with uniform density over the top surface of a low ther- 
mal conductivity material whose bottom surface is in good 
thermal contact with a heat sink. The only path for heat flow is 
through the low thermal conductivity material. The function 
^(;ic) is the height of the top surface of that material, and the 
temperature profile at that surface is proportional to y{x). 

Figure 3 illustrates the different effects of lateral 
thermal gradients on the thermal contrast displayed 
by the conventional design and the new design. 
For this comparison, it is assumed that 1) the same 
lateral gradient exists in the temperature T^(x) of 
the back plate of the conventional target and in the 
temperature Ts(x) of the substrate of the new type 
target, 2) the temperature of the back plate of the 
conventional target is controlled to exceed the tem- 
perature of the front plate by a fixed difference at 
the point x =0, and 3) the power dissipated in the 
new type target is set to give that same temperature 
difference at the points of discontinuity in the 
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height profile y(_x). The striking feature is that the 
distortion of the temperature profile 7'(jc) gener- 
ated by the conventional target deviates from the 
ideal much more than that of the temperature pro- 
file T[y(x)] generated by the new type of target. 

The next section of this paper presents the the- 
ory of operation of a thermal test target based on 
the new design. That section is followed by one 
that describes some targets built to demonstrate the 
new design, a section that compares the measured 
performance of these devices with the theoretical 
predictions, and a final section that illustrates how 
the targets would be used to extend the contrast 
range of a conventional target. 




Figure 3. Comparison of effect of temperature gradients on new 
and conventional approaches to thermal resolution targets. Part 
a) applies to a conventional target. Tr(x) is the temperature of 
the front plate between the slots, Tb(x) is the temperature of the 
back plate, and T(x) is the temperature profile generated by the 
target. Part b) applies to a target based on the design described 
in this paper. The liinction y(x) is the height profile of the top 
surface of the target, Ts(x) is the temperature gradient in the 
target substrate, and T[y(x)] is the temperature profile gener- 
ated by the target. 

2. Theory of Operation 
2.1 The Model 

An idealized device is used to model the heat 
transfer and temperature distribution; its cross sec- 
tion simulates a thin slice of material of the actual 



device which is shown in view B-B of figure 5. The 
symmetry properties of this device were used ex- 
tensively in this model. Figure 5 shows the mirror 
symmetry about each center line. Consequently, 
figure 4 extends only from the center line to the 
edge along section-line B-B. Because mirror sym- 
metry exists about the center line B-B, a tempera- 
ture profile for this cross section is obtained using a 
two-dimensional model. 
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Figure 4. An idealized device for computer modeling of various 
sources of error affecting the new type of thermal resolution 
target described in this paper. 
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Figure 5. Schematic diagram of the fabricated target (not drawn 
to scale). All dimensions are in millimeters. The views A-A and 
B-B are cross sections of thin slices taken at the positions indi- 
cated. The heights of these cross-sectional views are exagger- 
ated to illustrate the inconel and aluminum fdms and the shape of 
the groove. The X and Y labels correspond to those in figure 4. 
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The following assumptions were made for this 
model: 

1) The bottom surface of the device is in contact 
with a perfect heat sink, and all temperatures 
are measured relative to the temperature of the 
heat sink; 

2) The top surface of the device has three heaters, 
a groove heater, a center heater, and an edge 
heater, and they can be independently heated; 

3) The sides and the top surfaces above the three 
heaters are composed of perfect insulators; 

4) The heat transfer from the top surface of the 
device to the heat sink is entirely by conduction 
through the device; and 

5) Steady-state conditions exist; the temperature is 
not a function of time. 

The model is a finite difference model and con- 
sists of the cross-sectional area of figure 4 being 
divided into cells. A sketch of this cell construction 
is shown in Appendix 1. All x-cell dimensions (the 
lengths of the cell boundaries) are equal, and all 
^'-cell dimensions (the heights of the cell 
boundaries) are equal. The temperature of each cell 
is that at its center. The heat flow equations are 
written in integral form for each family of cells 
described below; 

Interior cells — The cell dimensions are xhy y. The 
net flow of heat into the cell is zero which means 
that the total heat flowing into the cell is equal to 
the total heat flowing out of the cell. As an exam- 
ple, the equation for the net heat flow into this cell 
per unit length, F{m, «), is given by 

0=F{m,n) = k{{x/y)[{T[m,n + \]-T[m,n]i) 
+{nm,n-\]-T[m, n])\+{y/x)\(T[m 
-\,n-\-T\m,n])+{T\m + \,n] 
-r[m, «])]}, 

where: 

T[m , « ] is the temperature of the m , n interior cell, 
k is the thermal conductivity of the material in 

watts per unit length per degree kelvin, 
X or ;; in the numerator is the length of the cell 

face across which heat flows, and 
X or ;; in the denominator is the distance between 

the centers of two adjacent cells. 

Surface cells with heaters — The cell dimensions are 
X by y/2. The heat flow into this cell from the 
heater, F(m, n), is equal to the heat flow across the 



lower horizontal cell boundary plus the net heat 
flow across the two vertical cell boundaries. Note 
that here F(m, «) is the heat per unit length 
supplied by the heater and not the net flow into the 
cell and is not zero. 

Exterior corner cell such as the cell of point (xi, y-i) in 
figure 4 — The cell dimensions are x/2 by y/2. The 
heater on the horizontal surface of this cell is half 
the size of those on the surface cells; consequently, 
the heat flow from the heater is only half that from 
the surface cell heater. This heat flow into the cell 
is equal to the heat flow across the lower boundary 
of the cell plus the heat flow across the vertical 
boundary of the cell. Because the heater and cell 
boundaries are both one-half of the dimensions of 
those for the interior cell, the heat flow equation 
obtained is the same as that for an interior cell. 
Again, F{m, «) is the heat per unit length supplied 
by the heater and not the net flow into the cell and 
is not zero. 

Interior corner cell such as the cell of point (x^, yj of 
figure 4 — The cell dimensions going clockwise 
around the cell, starting with the heater dimension 
are x/2 by y/2 by x/2hyyhy x by y/2. The heat 
flowing into this cell from the heater, which is half 
the heat flowing into a surface cell, is equal to the 
net heat flowing across the vertical boundaries plus 
the net heat flowing across the horizontal 
boundaries. The heat flowing across the surface 
represented by the first >'/2 dimension of this cell is 
zero. 

Boundaries across which no heat flows because 
of the constraints imposed by this model (assump- 
tion 3 above when no heater is present on the sur- 
face) are handled by setting the average 
temperature gradient across this surface equal to 
zero. This is done, for the algorithm given in the 
program in Appendix 1, by setting up virtual cells 
adjacent to the surface cells and then setting their 
temperatures the same as those of the adjacent 
(now interior) surface cells. 

Surface cells with heater elements are handled 
by this algorithm by setting up virtual cells as be- 
fore. The heat dissipated by the heaters is doubled 
by the program to accommodate the imaginary 
heat flow into the virtual cells and thereby provide 
the specified heat flow into the surface cells under 
the heaters. 

The equations obtained by applying all of the 
cell conditions stated are solved for the tempera- 
ture of the cell. The program source code is listed 



634 



Volume 95, Number 6, November-December 1990 

Journal of Research of the National Institute of Standards and Technology 



in Appendix 1 and can be compiled by version 3 of 
TURBO PASCAL.'' ^ The program first sets all the 
temperatures to zero, the temperature of the heat 
sink. It then applies the heat dissipation specified 
for the heaters and successively calculates the tem- 
peratures of all cells. This process is repeated 500 
times. Every successive iteration uses the tempera- 
tures calculated in the previous iteration with the 
same specified heater power dissipation. The tem- 
peratures of all the cells converge to within a small 
fraction of 1% of their final values in this many 
iterations. 

This program also allows the width and the 
depth of the groove to be independently varied in 
discrete steps, the power density in the edge heater 
to be set to zero or to the same value as that in the 
center heater, and the power density in the groove 
heater to be set to be a fraction between zero and 
the value set for the center heater. A typical result 
with the power density in the groove and edge 
heaters set the same as that in the center heater is 
shown in figure 6. The power density per unit con- 
ductivity in the center heater was adjusted to nor- 
malize the temperature at the top surface to 100 for 
these calculations. The same value was used in 
later calculations having the same geometry. 
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Figure 6. A typical result of the computer program listed in 
Appendix 1 for the device of figure 4 when the groove and edge 
heaters are dissipating the same power density as the central 
heater. For this example, there are 25 nodes on the x-axis, 11 
nodes on the ^i-axis, >'i=0.5 mm, ^2 = 1-0 mm, Xi=2.5 mm, 
X2=3.33 mm, and X3 = 5 mm, the lower left node is at the point 
(0 mm, mm), and the upper right node is at the point (5 mm, 
1 mm). 



' Certain commercial equipment, instruments or materials are 
identified in this paper to specify adequately the experimental 
procedure. Such identification does not imply recommendation 
or endorsement by the National Institute of Standards and 
Teachnology, nor does it imply that the materials or equipment 
identified are necessarily the best available for the purpose. 
^ TURBO PASCAL is a registered trademark of Borland Inter- 
national, Scotts Valley, CA. 



The important point of figure 6 is that the tem- 
perature profile of the top surface follows the phys- 
ical profile of the top surface. This result is a direct 
consequence of the assumptions used in deriving 
the model, and of the assumption of uniform dissi- 
pation of the same power density in all heaters. 
How well these assumptions can be realized in 
practice will, in part, determine how well the per- 
formance of a real device can approach the perfor- 
mance of the idealized device. The next three 
figures show some of the deviations from idealized 
performance that can be expected with different 
power densities in the device heaters. 

Figure 7 shows the effect of groove depth on the 
temperature profile (with the same power density 
dissipated in the center and edge heaters and no 
power dissipated in the groove heater) for a target 
having the geometry shown in figure 4 and the 
same groove width and overall width as the target 
analyzed in figure 6. This figure shows that the 
temperature profile no longer follows the physical 
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Figure 7. Normalized surface temperature for the device of fig- 
ure 4 when the edge and center heaters are dissipating the same 
power density, but the groove heater is dissipating no power. 
The dimensions are the same as those of the device modeled in 
figure 6 except that the cases of yi=0.2 mm (solid line), 
^'1=0.5 nun (dashed line), and ;'i=0.8 mm (dot-dash line) are 
considered. 
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profile of the top surface of the target, and that the 
largest deviation fi-om uniformity on each side of 
the temperature discontinuity occurs with shallow 
grooves. Therefore, for sharp thermal contrast, the 
groove heater must dissipate about the same power 
density as the center heater, and the target must 
have grooves that extend almost to the bottom of 
the low thermal conductivity material. The appli- 
cations of these criteria are depicted in figure 8 
which shows that fairly good uniformity on each 
side of the temperature discontinuity can be 
achieved with a groove depth that is half of the 
thickness of the low conductivity material and a 
power density in the groove that is half of that in 
the center heater. 
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Figure 8. Normalized surface temperature for the device of fig- 
ure 4 with the dimensions of the device modeled in figure 6 
when the groove heater is dissipating half the power density 
being dissipated by the center and edge heaters. 

Figure 9 shows the effect of dissipating no heat 
in the edge heater while maintaining the same 
power density in the groove and center heaters. 
The loss of uniformity is striking, but a sharp dis- 
continuity is still evident, even though the magni- 
tude of the discontinuity is only 60% of its value in 
figure 6. The poor uniformity shows that, with no 
edge heater, the center heater should cover the 



entire top surface of the low conductivity material 
in a target of this type. 

In the next section, devices are described that 
were built to test the criteria developed in this 
section. 
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Figure 9. Normalized surface temperature for the device of fig- 
ure 4 with the dimensions of the device modeled in figure 6 
when the groove and center heaters are dissipating the same 
uniform power density, but the edge heater is off. 

3. Device Fabrication and Testing 

A top view and two cross sections of one of the 
targets fabricated to test these criteria are shown in 
figure 5. In this target, there are no edge heaters, 
and the groove and center heaters are connected to 
the same aluminum contact pads. 

Fused silica (amorphous Si02) was used as the 
low conductivity material because its thermal con- 
ductivity is low compared to that of aluminum, but 
high compared to that of air. A 10-nm-thick layer 
of sputter-deposited inconel was used as the surface 
heater because inconel has a relatively high resis- 
tivity, and this thickness of inconel film produced a 
heater with both a usable resistance between 100 
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and 250 O and a reasonably high mfrared emissiv- 
ity. Figure 10 shows the spectral dependence of the 
measured reflectance and transmittance, and the 
calculated emissivity, of a nominal 10-nm film of 
inconel on a 750-/xm- (0.030-in-) thick piece of 
fused silica. Aluminum was used for the heater 
contact pads. The fused silica device was attached 
to the heat sink with a fast-curing epoxy adhesive. 
Aluminum was used as the heat sink material be- 
cause it is an inexpensive material that is easily ma- 
chined and has a high thermal conductivity. 
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Figure 10. The emissivity e (dashed line) of the inconel-fused 
silica structure used in fabricating the devices of figure 5. The 
emissivity was calculated from e = \—r—t, from measured data 
for the reflectance r and the transmittance t of the structure. 
These data are also shown in the figure. 

Table 1 compares the thermal conductivities of 
aluminum, fused silica, and air. Table 2 compares 
the electrical resistance per square of a lO-nm film 
of inconel and a 250-nm film of aluminum. Finally, 
notice that the net power density per unit tempera- 
ture difference radiated by a black-body at 300 K is 
only about 0.5% of that conducted through a 1-mm 
thick piece of fused silica. All of these results sug- 
gest that the assumptions stated in connection with 
figure 5 are well approximated by this device. 

The fused silica was cut and polished using con- 
ventional optical shop techniques into 75-mm- (3- 
in-) diameter, 750-;xm- (0.030-in-) thick wafers for 
subsequent processing. Rectangular grooves whose 
long axes were parallel to the future current direc- 
tions were ultrasonically machined into the top sur- 
face of the fused silica wafers. The walls on the 
long axis of the groove were vertical to produce an 
abrupt change in temperature across them (as 
shown in figs. 4 and 5); the walls at the ends of the 
grooves were sloped to allow continuous metal 
film coverage permitting the current to go down 
into and up out of the groove. 
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Table 2. Electrical resistance of fUms 


Material 


Nominal thickness 
(nm) 


Resistance 
(fl/O) 
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inconel 


250 nm 
10 nm 


0.14 
100. 



Photolithographic techniques were used to pro- 
duce the thin film heater and contact pads on the 
top surface of the fused silica. A number of differ- 
ent variations of the basic technique were tried in 
an attempt to find a technique that would allow the 
devices to be made without any special post-pro- 
cessing steps. Some approaches were found to be 
better than others, but no completely satisfactory 
approach was found. The presence of the grooves 
prevented the photoresist from behaving normally 
during spinning, curing, and etching. The ap- 
proaches tried and the results obtained are de- 
scribed in Appendix 2. 

After processing, the devices were cut from the 
wafers using standard techniques for fused silica. 
(Attempts to cut the wafers with a saw designed 
for dicing silicon wafers resulted in a great deal of 
wafer cracking, and breakage further limited the 
number of devices for testing. No such problems 
were encountered with the standard techniques for 
fused silica.) After cutting, the devices were 
mounted with a thin film of epoxy on 5-cm by 5-cm 
by 3-mm (2-in by 2-in by 1/8-in) heat sinks 
machined from soft (for higher conductivity) 
aluminum. 

Two bare copper wires were used as heater 
leads. Two layers of epoxy were used to attach 
each wire near one end of the heat sink to provide 
strain relief while insulating the wire from the heat 
sink. Then one end of each wire was attached to a 
device contact with silver-filled lacquer. These fab- 
ricated devices were then tested for electrical 
continuity. 
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4. Device Performance 

A PtSi camera operating in the 3- to 5-;im ther- 
mal infrared was used to study the performance of 
one of the low-contrast thermal resolution targets 
described above. The output of the camera was 
available on a television monitor as a thermal im- 
age and on an oscilloscope as a voltage vs time 
graph of one of the individual scan lines making up 
the television image. 

First, the camera was focused on a portion of a 
conventional 10-cm-square (4-in-square), four-bar 
thermal resolution target. The thermal contrast be- 
tween the rear heater and the slotted sheet was ad- 
justed to zero, and the image was recorded and 
stored for background subtraction. The contrast of 
the target slots was then adjusted to 4 K, the image 
recorded, the background subtracted, and the gain 
of the camera-oscilloscope system calibrated. Fig- 
ure 11, which does not reproduce well, shows the 
image of the target presented on the television 
screen during this calibration. Figure 12 shows an 
oscilloscope trace of one of the scan lines in the 
central part of the image shown in figure 11. 



The camera was then directed toward the low- 
contrast target, and the latter was observed with 
zero voltage across the heater to record the back- 
ground image. The background was observed to be 
unstable. The cause of the instability was traced to 
an image of the camera operator being inadver- 
tently reflected into the camera's field of view by 
the nonzero reflectance of the target. This points 
out the desirability of coating the top surface of the 
target with a thin, low-reflectance (high-emissiv- 
ity) material such as gold black [3], not only to de- 
crease the power required for a given contrast, but 
also to reduce the sensitivity to the background. To 
solve this problem for the existing devices, a black 
cloth was suspended around the target and camera, 
and a stable background was obtained. The latter 
was then recorded and stored for background sub- 
traction. 




Figure 11. Television image of a portion of the conventional four-bar target used for calibrating the gain of the PtSi CCD- 
array camera. The bright lines are the television-scan lines. The thermal contrast between the bars and the background was 4 K, and 
the background was subtracted to give the dark field background shown. 
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Figure 12. A scope trace of one of the scan lines in the image 
shown in figure 11. Nominal averages of the high- and low-tem- 
perature portions of the target are shown by the broken lines. 
The difference between the lines is about 75 IRE units. 



The voltage across the heater was then increased 
until the contrast between the heated top surface of 
the target and the surrounding heat sink was about 
4 K as measured on the calibrated oscilloscope, and 
the heater power was recorded. Figure 13, which 
also does not reproduce well, shows the image ob- 
tained under this condition. Figure 14 shows an os- 
cilloscope trace of one of the TV raster-scan lines 
of the camera. This particular line is located near 
the central part of the image shown in figure 13. 



heater. The rounding at the edges of the groove is 
partially an artifact caused by not assuring that the 
image of the groove exactly filled an integral num- 
ber of pixels on the PtSi array in the camera. This 
was verified to be the case by translating the target 
perpendicular to the optical axis of the camera with 
a micrometer-driven translation stage. It was possi- 
ble to make the pixel on either side of the groove 
assume any value between that in the center of the 
groove and the edge of the groove as the target 
was translated. 




Figure 14. A scope trace of one of the scan lines in the image 
shown in figure 13. Nominal averages of the high- and low-tem- 
perature portions of the target are shown by the broken lines. 
The difference between the lines is about 13 IRE units. 




Figure 13. Television image of the new type of low-contrast 
thermal resolution test target obtained in the same way as that 
shown in figure 1 1 . 

The rounding of the temperature profile in figure 
14 at the outside edges of the target was expected, 
as shown in the last section, because there was no 



The cause of the temperature gradient across the 
target is not known. A variation in the thickness of 
either the epoxy between the heat sink and the 
target or of the inconel heater film on the top sur- 
face of the target could cause this. Since no precau- 
tions were taken to assure the uniformity of the 
epoxy joint between the target and heat sink, a 
variation in thermal resistance caused by a varia- 
tion in epoxy thickness is the most likely cause for 
the temperature gradient observed with the cam- 
era. Since the thermal conductivity of the epoxy is 
about one-third that of fused silica, the thermal re- 
sistance of the 0.075-mm-thick layer of epoxy is 
about 33% of that of the 0.75-mm-thick fused silica 
layer. Therefore, a thickness gradient of the order 
of half the nominal thickness of the film would be 
needed to explain the observed temperature gradi- 
ent. This is not unreasonably large considering that 
no attempt was made to obtain a uniform joint. 
Clearly, it would be desirable to devise a way to 
assure a uniform epoxy joint. 

Figure 15 compares the theoretical and experi- 
mental temperature profiles for the device shown 
in figure 5. The experimental profile was obtained 
by reflecting the left-hand side of the oscilloscope 
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trace of figure 14 about the center of the groove 
and averaging it with the right-hand side of the 
trace. The theoretical profile was obtained using 
the model described earlier. Since this model does 
not account for the epoxy joint between the target 
and the heat sink, the device was modeled as hav- 
ing a thickness of 1.00 mm instead of the measured 
thickness of 0.75 mm. This additional thickness 
gives the same thermal resistance between the top 
surface of the device and the heat sink as 0.75 mm 
of fused silica in thermal series with 0.075 mm of 
epoxy. 
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Figure 15. Comparison of two parameters fit (filled circles) of 
theoretical model to the experimental data (open circles) shown 
in figure 14. The data for the left- and right-hand sides of target 
were averaged to obtain the experimental data plotted here. 

The model was fitted to the measured data by 
adjusting the groove depth until a good fit was ob- 
tained. The fit shown in figure 15 was obtained 
with a model groove depth of 0.20 mm, whereas 
the actual depth was 0.30 mm. The deeper groove 
on the actual device indicates that the actual con- 
trast obtained is not as good as that predicted by 



the model. This could be associated with a nonuni- 
form epoxy film or a small lateral thermal gradient 
existing in the heat sink under the groove. It might 
also be associated with a heat loss from the vertical 
groove wall of the actual device; a heat loss of this 
tj'pe would tend to reduce the contrast. Despite 
this unresolved problem with the fit of the calcu- 
lated temperature profile from this very simple, 
ideal model to the data from the actual devices that 
we built, figure 15 makes it clear that all of the 
major effects are explained by this model. Other 
than unwanted temperature gradients, there should 
be little or no deviations from the model in opti- 
mized targets. The target that we designed was 
built before the model calculations were devel- 
oped. Consequently, the target was used to test the 
model and provide insight for future target designs. 
Simple improvements in design (such as a groove 
depth more nearly equal to the target thickness and 
a heater covering the entire top surface of the 
target) should yield much more ideal devices. 

5. Use of New Type Target with 
Conventional Target 

The lower limit of thermal contrast available 
from a conventional target can be extended to ther- 
mal contrasts as low as those that can be resolved 
by any given camera with one of the new types of 
targets and an oscilloscope. The procedure is the 
following: 1) the camera/oscilloscope system is 
calibrated with the conventional target generating 
a reliably large thermal contrast, 2) the calibrated 
camera/oscilloscope system is used to calibrate the 
new type of target at the same order of thermal 
contrast, 3) the power in the new type of target is 
reduced to produce the desired thermal contrast, 4) 
the reduced thermal contrast is calculated using the 
formula presented below, and 5) the gain of the 
camera/oscilloscope system is increased and cali- 
brated. 

The formula for this type of calibration is 






(1) 



where DT^ov, is the lower level of thermal contrast 
generated in the new type target, i'low is the power 
dissipated in the new type target to generate DT\owy 
Z)i?cai is the difference between the in-groove and 
out-of-groove output from the camera/oscillo- 
scope when the power Pcai> i'low is dissipated in the 
new type target, DT^ is the thermal contrast from 
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the conventional target, and DR^ia is the difference 
between the in-slot and out-of-slot output from the 
camera/oscilloscope when the camera is pointed at 
the conventional target. It is important to notice 
that the gain of the camera/oscilloscope system 
must be the same for the traces from which DR^ 
and DRcai are read. Only after the calibration of the 
new type target at the higher power level can the 
camera gain be changed. 

Figures 12 and 14 illustrate how the target 
shown in figure 5 could be used as a low-contrast 
target. The thermal contrast of the conventional 
target was set to 4 K for the camera/oscilloscope 
trace shown in figure 12. Therefore, DTsui/DRm is 
about 4 K per 75 IRE units' =0.053 K per IRE 
unit. For the oscilloscope trace shown in figure 14, 
the voltage drop across and the current passing 
through the target heater were 8 V and 44.5 mA, 
respectively, and the thermal contrast between the 
groove and trace outside the groove is about 13 
IRE units. Therefore, DR^/P^= 15 IRE units per 
(8 VX44.5 mA)=42 IRE units per W. After the 
gain of the camera/oscilloscope system was in- 
creased and the voltage across the heater of the 
new type target decreased until the groove was 
just detectable in the image of the target, the 
voltage across and the current passing through the 
target heater were 3 V and 17.5 mA, respectively. 
(It is interesting that the groove could not be de- 
tected in the oscilloscope trace at this level, even 
though it could be detected in the thermal image.) 
Therefore, the minimum resolvable temperature 
difference is P,ow=52.5 mW, and DTi„^=(0.053 K 
per IRE unit) X (42 IRE units per watt) X 0.0525 
W=0.12 K. 

This minimum resolvable temperature difference 
of 120 mK, measurable with the PtSi camera re- 
ported above, is in good agreement with the value 
determined through very careful measurements 
against the conventional thermal resolution target. 
Why then is there a need for a new type of low- 
contrast target? The need is for the calibration of 
future cameras with better than 120-mK resolution. 
These experiments show that the new target is ca- 
pable of calibrations to near K. This is because 
the noise decreases linearly with the contrast or 
signal, and the signal-to-noise ratio remains con- 
stant allowing contrast calibrations to near K. 
This would be impossible, or at least very difficult 
to do, with the conventional target. 



^ The vertical deflection of the oscilloscope trace, as seen in 
figures 12 and 14, are given in IRE units. 
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6. Appendix 1. Source Code for Temperature Distribution Program 



PROGRAM compute_lctrtt_temperature_distribution; 

{USES Crt;) {REMOVE THE CURLY BRACKETS AROUND "USES Crt;" TO COMPILE IN 

TURBO PASCAL 4 AND TURBO PASCAL 5 . ) 
VAR T : ARRAY[-1. .26,0. .11] OF REAL; 

n, m, 1, heater_edge , Ml, M2, M3 , SUCC_M3, Nl, N2 , SUCC_N2 : INTEGER; 
answer : STRING [1] ; 
X, y, flux_per_k, fraction : REAL; 
{ < groove I central heater | edge | 
heater heater 

N2 *-4--*-4--*-4--*-4--*-4--*-4--*-4--* 

■+ * point (m,n) 



+ - 
* 

+- 



Nl *-4--*-4--*-4--* 

* I * I * I * 
I---I---I 



* 
* 



* 
* 
* 



* 
* 
* 



* 
* 
* 



* 
* 



* 

- + 

* 

■ + 
* 

■ + 



I---I 
I I 



cell 



*-4-_*_4-_*_4--*-4--*-4--*-4--*_4--*_4-_*.4-.*_^ — * 

Ml M2 M3 

I X I 



-+-*-+-* surfaces 

+ 
* 

+ 



) 



PROCEDURE initialize; 

BEGIN 

CLRSCR; 

WRITELN; 

WRITELN; 

WRITELN( 

' The lower left cell is indexed as (0,0) in this finite cell model.'); 

WRITELNC There are three heaters on the top surface: 1) groove heater, '); 

WRITELNC 2) central heater, ') 

WRITELNC 3) edge heater. '); 

WRITELN; 

WRITEC Index < 26 of x-axis cell at right edge of device: '); READLN(M3) ; 

SUCC_M3 := SUCC(M3) ; 

WRITEC Index of x-axis cell at edge of central heater: '); READLN(M2) ; 

WRITEC Index of x-axis cell at edge of groove: '); READLN(Ml) ; 

WRITEC Index < 11 of y-axis cell at top surface of device: '); READLN(N2) ; 

SUCC_N2 := SUCC(N2) ; 

WRITEC Index of y-axis cell at top of groove: '); READLN(Nl) ; 

FOR m := -1 TO SUCC_M3 

DO FOR n := TO SUCC_N2 

DO T[m,n] := 0; 
WRITEC Width X of device in mm: '); READLN(x) ; x := x/M3 ; 
WRITEC Height Y of device in mm: '); READLN(y) ; y := y/N2; 
WRITEC Power density per conductivity in central heater: '); 
READLN(flux_per_k) ; 
WRITEC Power density in groove heater relative to that in central heater: 



'); 
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READLN( fraction) ; 

WRITEC CR turns edge heater on, any other key turns it off: '); 
READLN( answer) ; 

IF answer = '' THEN heater_edge := M3 
ELSE heater_edge := M2; 
END; 

PROCEDURE list_temperatures; 
BEGIN 
CLRSCR; 

FOR n := N2 DOWNTO 
DO FOR m := TO M3 
DO BEGIN 

IF ( (ni>=Ml) OR (n<=Nl) ) THEN WRITE(T[m,n] :3 :0) 

ELSE WRITEC ') ; 
IF m - MS 
THEN BEGIN 

WRITELN; 
WRITELN; 
END; 
END; 
END; 

PROCEDURE adjust_temperatures; 

VAR XX, yy, bottom, driver : REAL; 

BEGIN 

XX := x*x; 

yy :- y*y; 

bottom := (xx + yy) ; 

driver := flux_per_k*x*y /bottom; 

FOR n := N2 DOWNTO 1 

DO BEGIN 

IF n > Nl 

THEN T[Ml-l,n] := T[Ml+l,n] 

ELSE T[-l,n] := T[l,n] ; 

T[M3+l,n] := T[M3-l,n] ; 

FOR m := TO MS 

DO BEGIN 

IF m < Ml 

THEN T[m,Nl+l] := T[m,Nl-l] 

ELSE T[m,N2+l] := T[m,N2-l] ; 

IF ( (n=Nl+l) AND (m=Ml) ) 

THEN T[Ml-l,n] :=T[Ml+l,n]; 

IF NOT ( ( (m<Ml) AND (n>Nl) ) OR ( (m=Ml) AND (n=Nl) ) ) 

THEN BEGIN 

T[m,n] := yy*(T[m-l,n]+T[m+l,n]) + xx*(T[m,n-l]+T[m,n+l] ) ; 
T[m,n] := T[m,n]/(2*bottom) ; 
IF n = N2 

THEN IF ( (m>=Ml) AND (m<=heater_edge) ) 
THEN T[m,n] := T[m,n] + driver 
ELSE 
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ELSE IF n - Nl 

THEN IF m <= Ml 

THEN T[m,n] := T[m,n] + f raction*driver ; 
END 
ELSE IF ( (m=Ml) AND (n=Nl) ) 
THEN BEGIN 

T[ra,n] := 

yy*T[m-l,n]+2*yy*T[m+l,n] + 2*xx*T[m,n-l]+xx*T[m,n+l] ; 
T[m,n] := T[m,n]/(3*bottoin) ; 
T[m,n] := T[m,n] + fraction*driver/3; 
END; 
END; 
END; 
END; 

BEGIN {Main program) 

initialize; 

FOR 1 :=-= 1 TO 500 

DO adjust_temperatures ; 

list_temperatures ; 

REPEAT UNTIL KEYPRESSED; 

END. 



7. Appendix 2. Fabrication of 
the Test Targets 

The targets were fabricated on 75-mm- (3-in-) 
diameter, 750-ju.m-tliick, fused silica wafers. In this 
fabrication process, grooves were first ultrasoni- 
cally machined into the wafers at specified loca- 
tions. Inconel-600 films'* 10 nm thick were sputter 
deposited on the wafers. Over the inconel films, 
pure aluminum films 250 nm thick were sputter-de- 
posited. The target patterns were aligned to the ul- 
trasonically machined grooves and patterned using 
photolithographic techniques. Only a small number 
of machined wafers were available which severely 
limited the amount of experimentation that could 
be done for process development. This number was 
further limited because the bottoms of the ultrason- 
ically machined grooves were rough enough to 
prevent good electrical continuity in the 10-nm- 
thick inconel film deposited in them. Attempts to 
smooth these machined surfaces were made by dip- 
ping the wafers into 49% HF. After a few minutes 
etching, sharply outlined surface scratches ap- 
peared on the silica surfaces rendering some wafers 
useless for target fabrication. 



The fabricated targets are shown schematically 
in the dimensioned drawing of figure 5. Side-by- 
side arrays of six of these targets were simulta- 
neously fabricated on each wafer. The outer two 
targets on each side of the array, i.e., the targets 
closest to the periphery of the wafer, contained no 
grooves and served as control targets during fabri- 
cation. 

Three different methods were used to fabricate 
these targets. One of these methods consisted of a 
lift-off technique which is described later. The 
other two consisted of differential etching which 
made use of the fact that the aluminum etchant 
etched inconel much slower than aluminum. This 
enabled the uppermost aluminum film to be pat- 
terned photolithographically and etched off of the 
inconel film. The primary differences in these two 
etching methods are shown in table A2.1. The re- 
sistances, in ohms, of the target heaters, fabricated 
by each method, are given in table A2.2. In this 
process, the first level was patterned into photore- 
sist and consisted of the aluminum busses at the 
ends of the target grooves (shown in fig. 9). The 
aluminum was etched at 45 °C with a commercial 
aluminum etchant.' The etching was stopped as 



* The inconel-600 has a composition of 75.6 wt.% Ni, 8.5% Fe, 
and 15.9% Cr. 



' This etchant consisted of 85 parts of 85% phosphoric acid, 2 
parts of 70% nitric acid, and 13 parts of glacial acetic acid. 
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soon as the aluminum film completely dissolved ex- 
posing the underlying inconel film. The resist was 
then stripped with acetone, and the second level, 
which consisted of the inconel heaters, was pat- 
terned into another photoresist film. This patterned 
photoresist film completely covered the aluminum 
busses protecting them fi-om attack by the inconel 
etchant. The inconel was etched with Cyantek CR- 
9 Chromium Photomask Etchant,* which is a con- 
centrated perchloric acid-based etchant. 

Table A2.1. Primary differences in the etching methods 



Differential etch method 1 



Differential etch method 2 



Sputter Deposition vacuum 
broken between the deposi- 
tions of inconel and aluminum 



Photoresist — Shipley Micro- 
posit 1350J 

Photoresist coated by conven- 
tional spinning and open hot 
plate baking 



Sputter Deposition with alu- 
minum deposited immediately 
over inconel without breaking 
vacuum 

Photoresist — Shipley Micro- 
posit 1470 

Photoresist coated by the 
flood-surface tension method 
and a "capped" hot plate bake 



Table A2.2. Resistances" of target heaters in ohms 







Target nimiber on 


Wafer 




Method 


l" 


2 


3 


4 


5 


6" 


Dif. etch method 1 


138 


33'' 


167 


150 


22-= 


136 


Dif. etch method 2 


163 


257 


244 


293" 


447d 


157 


Lift off 


183 


180 


197 


163 


171 


164 



' These resistances are measured by touching the probes of an 
ohmmeter to the aluminum busses. The resistances of the 
busses and probe contacts are assumed to be negligible. 

'' Have no grooves and serve as control targets. 

■^ These low values are attributed to the small strips of aluminum 
that were observed on the wafer surfaces along the grooves. 

■'Targets where the etching of the Shipley 314 developer was 
most severe. 



Difficulty was encountered with photoresist 
building up on the sharp inside (the sides toward 
the wafer center) precipices of the grooves when it 
was spun on the wafer by a conventional method. 
The photoresist in this buildup could not be satis- 
factorily patterned. A "flood-surface-tension" 
method was developed for spin coating the pho- 
toresist that did not produce this unwanted pho- 
toresist buildup. In this method: 1) the wafer is 
mounted on the spinner chuck and completely 



flooded with photoresist; 2) immediately following 
the flooding, the wafer is spun at 300 to 400 rpm 
for approximately 1.5 s, and then ramped to 1500 
rpm and spun for 5 s; 3) when the spinner stops, the 
wafer is allowed to set on the chuck for an addi- 
tional 60 s; and then 4) prebaked, in a flat, horizon- 
tal position, on a hot plate at 95 °C for 3.5 min with 
an aluminum foil tent covering it and the immedi- 
ate surrounding hot plate area. The purpose of the 
two short, slow spins is to remove excess photore- 
sist and achieve a reasonably thin resist coating 
without producing either buildup in or significant 
solvent loss from the photoresist. The 60-s setting 
period allows the photoresist surface tension to pla- 
narize the resist and reduce any buildup. The pre- 
bake is done with the wafer flat to prevent uneven 
resist film formation, and the alumimmi tent pro- 
motes better heating of the resist and improves its 
adhesion. It should be noted that the thermal con- 
ductivity of a 750-/xm- (30-mil-) thick fused silica 
wafer is substantially less than that of a silicon 
wafer and necessitated the use of the tent. 

Exposures of uv were made at 180 to 220 mJ/ 
cm^ using AR-chrome masks. Spray-puddle-devel- 
oping was used with full-strength, Shipley 314 
developer.^ Before exposures were made, the 
coated wafers were either allowed to set in the air 
for at least 1 h, or if the exposure had to be made 
sooner, the wafer was dipped in deionized water 
for 1 min and spun dry. The purpose of this step 
was to replace water in the photoresist film that 
was removed during the baking. This water is nec- 
essary for the resist photolysis to proceed. 

Some problems were encountered with each of 
the differential etch methods. With method 1, thin 
lines of unetched aluminum were observed on the 
wafer surfaces adjacent to the sides of the grooves 
toward the wafer center. These lines are believed 
to be a result of resist buildup. With method 2, pho- 
toresist was observed on the groove walls after de- 
veloping. The removal of this resist was attempted 
by re-exposing the groove walls and developing by 
puddling a drop of Shipley 314 developer in the 
grooves. The developer etched the aluminum 
where it was puddled which, in turn, caused the 
inconel to be etched thinner in these areas. 

7.1 Lift-Off Method 

Inconel was first sputter-deposited onto the 
wafer which was then conventionally spin-coated 
with Shipley 1350J and prebaked 3.5 min at 95 °C 
on a hot plate. A uv exposure of 220 mJ/cm^ was 



' A product of Cyantek Chemicals, Mountain View, CA. 



' A product of Shipley Co., Inc., Newton, MA. 
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used to pattern a mask of the aluminum busses (the 
polarity of which was the reverse of the masks 
used in the differential etching methods) onto the 
wafer. After this exposure, the wafer was soaked in 
chlorobenzene for 5 min and spray-puddle-devel- 
oped. A 250-nm aluminum film was evaporated, in 
a filament evaporator, onto the patterned wafer un- 
der a vacuum where /'<4xl0~' Torr. The dis- 
tance between the filament and the wafer was 
about 50 cm. The filament was a trough filament 
into which a few pieces of braided tungsten wire 
were placed with aluminum pellets. When the alu- 
minum melted, it wetted the tungsten wire which 
then gave even evaporation of the aluminum for 
the 250-nm-thick film. Lift-off of the aluminum film 
was then accohiplished by soaking the wafer for 3 
h in acetone with intermediate acetone flushes 
from a squeeze bottle to help dislodge the film. The 
aluminum in the corners of the grooves at the base 
of their vertical walls required "swabbing" with a 
sharp plastic point und/jr acetone for complete re- 
moval. 

The remaining inconel film was then patterned 
by coating the wafer with Shipley 1470 photoresist 
using the flood-surface-tension method, patterning 
the inconel mask with a 220 mJ/cm^ exposure, 
spray-puddle-developing with Shipley 314 devel- 
oper, and etched with CR-9. The resistances of 
these inconel heaters are given in table A2.2. 
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